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Abstract 

The General AntiParticle Spectrometer (GAPS) is a novel approach for indirect dark matter searches that exploits 
cosmic antideuterons. GAPS utilizes a distinctive detection method using atomic X-rays and charged particles from 
the exotic atom as well as the timing, stopping range and dE/dX energy deposit of the incoming particle, which pro- 
vides excellent antideuteron identification. Prior to the future balloon experiment, an accelerator test was conducted 
in 2004 and 2005 at KEK, Japan, in order to precisely measure the X-ray yields of antiprotonic exotic atoms formed 
with different target materials [I]. The X-ray yields of the exotic atoms with Al and S targets were obtained as ~ 75%, 
which are higher than were previously assumed in O. A simple, but comprehensive cascade model has been devel- 
oped not only to evaluate the measurement results but also to predict the X-ray yields of the exotic atoms formed with 
any materials in the GAPS instrument. The cascade model is extendable to any kind of exotic atom (any negatively 
charged cascading particles with any target materials), and it was compared and validated with other experimental data 
and other cascade models for muonic and antiprotonic exotic atoms. The X-ray yields of the antideuteronic exotic 
atoms are predicted with a simple cascade model and the sensitivity for the GAPS antideuteron search is estimated as 
1.3 xl0~ 6 [m _2 s _1 sr _1 (GeV/n) _1 ] for the proposed long duration balloon program [3 1, which suggests that GAPS has 
a strong potential to detect antideuterons as a dark matter signature. A GAPS prototype flight (pGAPS) was launched 
successfully from the JAXA/ISAS balloon facility in Hokkaido, Japan in summer 2012 JUS] and a proposed GAPS 
science flight is to fly from Antarctica in the austral summer of 2017-2018. 
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1. Introduction 

1.1. Overview 

An accelerator test was conducted in 2004 and 2005 
at KEK, Japan, in order to precisely measure the X-ray 
yields of antiprotonic exotic atoms formed with differ- 
ent target materials HI. This paper describes not only 
the detailed analysis of the X-ray yields for Al and S 
targets (Section B), but also the development of a com- 
prehensive cascade model for the exotic atom (Section 
|2j. It is crucial to develop the cascade model to evaluate 
X-ray yields of exotic atoms with a variety of target ma- 
terials since the GAPS project, an indirect dark matter 
search, utilizes atomic X-rays of exotic atoms to iden- 
tify antideuterons, and the exotic atoms can be formed 



with any materials in the GAPS instrument (see Sec- 
tion |1.4| >. As described in Section[2j the cascade model 
was compared and validated with the experimental data 
and other cascade models for muonic and antiprotonic 
exotic atoms. In Section |3.8| the cascade model was 
used to estimate the X-ray yields of antideuteronic ex- 
otic atoms formed in the GAPS instrument to derive the 
ultimate sensitivity for the antideuteron search. Note 
that even if antiprotons did not stop in the Si target, they 
could be identified with the atomic X-rays of the ex- 
otic atom formed with the known material in the GAPS 
instrument. The sensitivity is estimated as 1.3 xl0~ 6 
[nrV^r-^GeV/n)- 1 ] for the 105 day, 3 flight long du- 
ration balloon program 1 3 1 and thus GAPS has a strong 
potential to detect antideuterons as a dark matter signa- 
ture (see Section[T3|. 
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1.2. Dark Matter Candidates 

The recent result by the Planck experiment |6| shows 
that 68% of our universe is composed of dark energy, 
and 27% is dark matter (~ 5% for baryonic matter). The 
nature and origin of these phenomena, however, are still 
unknown, and thus are the great cosmological problems 
of the 21st century. Unlike dark energy, dark matter is 
well-motivated by many theoretical models, and many 
experiments are currently being conducted to determine 
the origin of dark matter. 

The existence of dark matter was postulated by Fritz 
Zwicky in 1933 from the observation of the rotational 
speed of galaxies. The recent observations of gravita- 
tional lensing in the Bullet Cluster (two colliding clus- 
ters of galaxies), also indicate the existence of dark mat- 
ter 0. 

Since dark matter has never been directly observed, 
it is considered to interact with the Standard Model 
particles only by the weak force and the gravitational 
force as seen in rotational curves and gravitational lens- 
ing. The small density fluctuations seen in the cos- 
mic microwave background (CMB) [8| and the large 
scale structure of the present universe indicate that dark 
matter should be a non-relativistic and massive parti- 
cle (called cold dark matter). Moreover, it should be 
stable on a cosmological time scale to be observed in 
the present universe. Weakly interacting massive parti- 
cles (WIMPs) are the theoretically best-motivated can- 
didates among the variety of dark matter candidates. 
Neutralinos, the lightest supersymmetric partner (LSP) 
in supersymmetric theories, and Kaluza-Klein particles 
(LKP) and right-handed neutrinos (LZP) in extra di- 
mension theories are examples of popular WIMP can- 
didates. 

1.3. Antideuterons for Dark Matter Search 

There are dozens of experiments designed to search 
for particles associated with various manifestations of 
WIMP dark matter categorized into three types, parti- 
cle collider, direct search, and indirect search. The di- 
rect and indirect searches will measure the relic WIMPs, 
while the particle collider will try to create WIMPs. The 
detection methods and the background models for each 
search are different, but also complementary, helping to 
illuminate the nature of dark matter. 

The General AntiParticle Spectrometer (GAPS) is a 
novel approach for indirect dark matter searches that 
exploits cosmic antideuterons. WIMPs annihilate with 
each other and produce a variety of standard model par- 
ticles including electrons, positrons, gamma rays, an- 
tiprotons and antideuterons. 
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Figure 1 : Antideuteron flux at the top of the atmosphere, compared 
with the BESS upper limit (5), and GAPS and AMS sensitivity 0. 
The sensitivity for the AMS 5 year flight was estimated, based on 1101 
4]. The blue dashed line (LZP), black dotted line (LSP), and green 
dot-dashed line (LKP) represent the primary antideuteron fluxes due 
to the dark matter annihilations [11']. The red solid line represents the 
secondary/tertiary flux due to the cosmic-ray interactions I12II13|[T41 . 

Antideuteron production in WIMP- WIMP annihila- 
tions was proposed by Donato et al., in 2000 ||T31 [TBI . 
The antideuteron flux due to WIMP- WIMP annihilation 
(called primary flux) can be estimated based on the dark 
matter density profile of the galaxy, the WIMP- WIMP 
annihilation channel, the hadronization and coalescence 
model, and the propagation model. The primary an- 
tideuteron flux due to the WIMP- WIMP annihilation 
is shown in Figure [TJ (solid purple line: LSP with m x 
~ 100 GeV, dashed green line: LKP with m x ~ 500 
GeV, dashed blue line: LZP with m x ~ 40 GeV) ifTTl . 
The relatively flat peak is located at E ~ 0.2 GeV/n. 
The antideuteron flux due to the cosmic -ray interac- 
tions with the interstellar medium (secondary/tertiary 
flux, red dashed line) is also shown in Figure [Jj 1 121IT31 
114*1 . Unlike primary antideuterons, collision kinemat- 
ics suppress the formation of low-energy secondary an- 
tideuterons. Moreover, the interaction rate is drastically 
decreased since the flux of the cosmic -ray protons fol- 
lows the power law, F p ~ E~ 2 1 . Therefore, the primary 
antideuteron flux is two orders of magnitude larger than 
the secondary/tertiary antideuteron flux at low energy, 
and we can clearly distinguish them. 

The GAPS and AMS (5 year flight) sensitivities 0, 
and the current upper limit for the antideuteron flux ob- 
tained by the BESS experiment [9| are also shown in 
Figure [JJ As seen in the figure, the GAPS experiment is 
more than two order of magnitude more sensitive than 
the BESS upper limit and 1.5 times more sensitive than 
the AMS satellite mission. (The sensitivity for a GAPS 
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210 day flight program (LDB+) is also shown in the 
figure.) Thus, GAPS has a strong potential to detect 
antideuterons as the dark matter signature. In the fol- 
lowing section, the details of the GAPS project are in- 
troduced including the detection concept and the instru- 
mental design. 

1.4. GAPS Project 

1.4.1. Overview of the GAPS Project 

The general antiparticle spectrometer (GAPS) was 
first proposed by Hailey and Mori in 2002 [2|, and was 
originally named the Gaseous AntiParticle Spectrome- 
ter. The original GAPS was designed to use a gaseous 
target, but with further studies, including the KEK (high 
energy accelerator research organization) beam test in 
Japan described below, we concluded that a solid tar- 
get was more efficient and effective for the flight ex- 
periment. GAPS is a balloon-borne experiment (flight 
altitude ~ 35 km), and there are constraints on the size 
and mass of the payload. Therefore, the solid target can 
greatly simplify the setup of the GAPS flight module by 
removing the bulky gas handling system and allowing 
more complex designs, such as a multi-layer tracker ge- 
ometry. The higher density of the solid target can also 
easily slow down and stop more incoming antiparticles, 
which provides a larger detectable energy range. 

A GAPS prototype flight (pGAPS) was launched 
successfully from the JAXA/ISAS balloon facility in 
Hokkaido, Japan in the summer of 2012 JU [5). The 
pGAPS flight instrument was composed of 3 layers of 
plastic scintillators and 3 layers of Si(Li) detectors (2 
detectors in each layer) inside the detector vessel be- 
tween the middle and the bottom plastic layers. Figure 
[2]shows the pGAS payload hung from the launcher (left 
picture), and the plastic scintillators and the Si(Li) de- 
tector vessel inside the gondola (right picture). About 
one million cosmic -ray events were triggered during the 
6 hours flight, (~ 3 hours at float, ~ 33 km). A proposed 
GAPS science flight is to fly from Antarctica in the aus- 
tral summer of 2017-2018. 

7.4.2. Detection Concept 

The GAPS detection method involves capturing an- 
tiparticles into a target material with the subsequent 
formation of an excited exotic atom. A time-of-flight 
(TOF) system measures the velocity (energy) and di- 
rection of an incoming antiparticle. It slows down by 
the dE/dX energy loss and stops in the target material, 
forming an excited exotic atom. The exotic atom de- 
excites in a complex process involving Auger ionization 
and electron refilling at high quantum number states, 




Figure 2: pGAPS payload in the launcher (left) and the plastic scintil- 
lators and the Si(Li) detector vessel inside the gondola (right). 

followed by the emission of X-rays at the lower quan- 
tum states (see Section [2]). With known atomic num- 
ber of the target, the Bohr formula for the X-ray en- 
ergy uniquely determines the mass of the captured an- 
tiparticle 0. Ultimately, the antiparticle is captured by 
the nucleus in the atom, where it is annihilated with 
the emission of pions and protons. The number of pi- 
ons and protons produced by the nuclear annihilation 
is approximately proportional to number of antinucle- 
ons, which provides an additional discriminant to iden- 
tify the incoming antiparticle. The entire process takes 
place in ~ 1 ns. This process is illustrated in Figure 
[3] The Lithium-drifted silicon (Si(Li)) detector will be 
used both as a target to form the exotic atom and as a 

/5 
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Figure 3: Schematic for GAPS detection method. An antiparticle 
slows down and stops in the Si(Li) target forming an exotic atom. The 
atomic X-rays will be emitted as it de-excites followed by the pion 
(and proton) emission in the nuclear annihilation. 
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detector to detect atomic X-rays and the annihilation 
products. The concept of the detection technique has 
been verified through the accelerator testing at KEK in 
2004 and 2005, as described in Section[3] 
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Figure 4: Antideuteron-antiproton identification technique: (1) depth 
sensing, (2) unique atomic X-rays, (3) pion multiplicity. 



Antiprotons are the major background in this exper- 
iment, since they can also form exotic atoms and pro- 
duce atomic X-rays and charged particles. However, 
the atomic X-rays and the number of pions and pro- 
tons emitted from the exotic atom uniquely identify the 
mass of the original antiparticle, as do the depth sensing 
(stopping range of the incoming particle) and the dE/dX 
energy loss in each Si(Li) detector, once the velocity of 
the incoming antiparticle is determined by the TOF sys- 
tem. 

The three highest antideuteronic X-rays with a Si tar- 
get in the GAPS detectable energy range are 67 keV, 
44 keV and 30 keV, while antiprotonic X-rays are 107 
keV, 58 keV, and 35 keV. The number of charged par- 
ticles produced by the nuclear annihilation for the an- 
tideuteronic exotic atom is approximately twice as many 
as the one for the antiprotonic exotic atom. Addition- 
ally, antideuterons with the same speed have a longer 
stopping range and can go deeper into the detector than 
antiprotons (see Figure|4|. Thus, antideuterons with the 
same stopping range will have a smaller velocity and de- 
posit more energy at each layer than antiprotons, since 
the dE/dX energy loss is inversely proportional to the 
velocity squared at low energy. As a result, these de- 
tection methods provide a greater than 10 6 antiproton 
rejection factor 131 . 



1.4.3. Instrumental Design 

The GAPS balloon flight will have several unique 
features. It will be the first balloon flight with a very 
large, pixellated Si(Li) detector surrounded by a very 
large TOF system without a pressure vessel (see Figure 
[5}. There will be 10 layers of detectors surrounded by 
TOF plastic scintillators, with each layer composed of 
4 inch diameter, 2.5 mm thick Si(Li) detectors. Each 
Si(Li) detector will be segmented into 4 strips, and ad- 
jacent tracking layers will have their strips positioned 
orthogonally, providing modest three-dimensional par- 
ticle tracking. 




Figure 5: GAPS detector design. 10 layers of Si(Li) detectors are 
surrounded by the TOF plastic scintillators. Each Si(Li) detector is 4 
inch diameter, 2.5 mm thick. 



A tracking geometry is a natural consequence of the 
current GAPS detection concept, unlike the original 
GAPS, whose cubic geometry only focused on detecting 
atomic X-rays using a gaseous target cell surrounded by 
scintillation detectors without determining trajectory in- 
formation Q][2][T7). The tracking geometry can count 
the number of particles produced in the nuclear annihi- 
lation and separately identify atomic X-rays from parti- 
cle tracks. It also permits direct measurement of parti- 
cle stopping depth and naturally conforms to the multi- 
detector geometry. Note that since each strip is rela- 
tively small, ~ 2 cm wide, X-rays and charged particles 
(pions/protons) can be detected separately in the differ- 
ent strips/channels. Each Si(Li) layer also works as a 
degrader and a target material to slow down the incom- 
ing antiparticle and to form an exotic atom. 
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2. Cascade Model for Exotic Atoms 

2.1. Overview of Cascade Model 

As seen in the previous section, X-ray yields of ex- 
otic atoms play an important role in the GAPS an- 
tideuteron detection. The energy of the atomic X-ray 
is unique to the exotic atom, allowing us to differen- 
tiate antideuterons from other particles, including an- 
tiprotons. Therefore, it is crucial to develop a compre- 
hensive cascade model to estimate the X-ray yields for 
any kind of exotic atom (any negatively charged cascad- 
ing particles with any target materials) that can form in 
the GAPS instrument. 

Cascade models for exotic atoms were widely devel- 
oped after the existence of the exotic atom was predicted 
in the 1940s. Since the GAPS project focuses on the an- 
tiprotonic and antideuteronic exotic atoms formed with 
variety of target materials, we have developed a gen- 
eralized and extendable cascade model. Additionally, 
since the GAPS detector is designed for X-rays with 
an energy higher than 10 keV, a very simple cascade 
model with a few parameters has been developed, fo- 
cusing on the low n state transitions (E >10 keV). The 
parameters were optimized by the measurement of an- 
tiprotonic exotic atoms with Al and S targets at KEK 
in Japan 2005. The extended cascade model was used 
to estimate the X-ray yields of the antiprotonic and an- 
tideuteronic exotic atom with a Si target and other ma- 
terials in the GAPS instrument to derive the ultimate 
antideuteron sensitivity (see Section[3]l. 

2.2. Cascade Transitions 

A negatively charged particle ([T , if , K~, p, d, etc., 
called "cascader" hereafter) will be captured into a tar- 
get atom at the radius of its outermost electrons after it 
slows down and its kinetic energy becomes comparable 
to the binding energy of an electron lfl8l [191 . The ini- 
tial principal quantum number for the exotic atom can 
be estimated as follows: 



n ~ n e ^M*/m* . 

Here, n e is the principal quantum number of the out- 
ermost electron shell of the target atom, m* is the re- 
duced mass of the electron in the target atom and M* is 
the reduced mass of the cascader. The cascade model 
is designed to calculate the probability for the cascader 
to be in the (n, I) state, where I is the orbital angular 
momentum, and to estimate the X-ray yields of the ex- 
otic atom as it decays. The cascade model starts at the 
electron K shell (n e = 1) and the orbital angular mo- 
mentum I is assumed to have a statistical distribution, 



Pi oc (21 + l)e al . There are (21 + 1) magnetic quantum 
numbers, m = — I + 1, -I + 2 ... ... I - 2,1 — 1, for each 
/, and e al is a correction factor due to the de-excitation 
at the outer shell, n e > 1 (a ~ 0.2 or less) |fl8l[T9l . The 
initial n in the cascade model is about 14 for [T, 16 for 
iC, 31 for K~, 42 for p, and 58 for d. 
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Figure 6: Schematic for the cascade model of the antiprotonic exotic 
atom. The Auger transitions dominate in high n states, while the ra- 
diative transitions dominate in low n states. The nuclear capture takes 
place in very low n states. 

The three leading de-excitation processes, Auger 
transition (emission of an Auger electron), radiative 
transition (emission of an atomic X-ray), and nuclear 
capture (interaction with the nucleus), dominate the cas- 
cade model for atoms with Z > 2, as shown in Figure [6] 
Auger transitions dominate at the beginning of the cas- 
cade, followed by radiative transitions. The nuclear cap- 
ture takes place in a very low n state. Since the exotic 
atom can be assumed to be hydrogen-like, the Auger 
and the radiative transitions with Al = ±1 dominate due 
to selection rules. 

2.2.7. Auger Transition 

In a high n state, an Auger electron is emitted as soon 
as the energy difference of the initial state (ri\, l\) and 
the final state (nz, h) exceeds the ionization energy. The 
Auger transition rate for the K shell and L shell elec- 
trons can be estimated by considering the interaction 
between the cascader and the electron as follows 



"l,'l->"2,'2 



32nac / Z* 
~Z 



a /i 

,.2 



max(Zi, I2) 



3(2/, + 1) 
y exp[y(4 arctan y - n)] ^ 
1 + y 2 sinh Try 
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pAKg.1 _ 16toc /Z*\ max(/i,/ 2 ) 

y 2 (4 + 5y 2 )(4 + 3y 2 ) 

(4+y 2 ) 3 
exp[y(4 arctany - n)] ^ 
sinh ny 

Here, p., y, and / are defined as follows. 



l> 

y 

T 



M/m e 



Z*a 



^(T/m e c 2 ) 2 + (2T/m e c 2 ) 
AE„, 



s f 

Jo 



drr , R(ni,h)R(n 2 ,h) 



j-Aug K r^Aug,L s j s ^ Auger transition rate for 
emitting K (L) shell electrons with the initial state (n\, 
l\) and the final state (/12, fe)> fl o is the Bohr radius of 
hydrogen atom, a is the fine structure constant, Z* is the 
effective nuclear charge seen from the electron, T is the 
kinetic energy of the emitted electron, and R(n, I) is the 
normalized radial function of the exotic atom. The tran- 
sitions with Al — +1 dominate the process, due to the 
transition selection rules as discussed above. Note that 
after the electrons are depleted by the Auger transition, 
the electrons can be filled from adjacent atoms with a 
refilling rate T re f and also from the higher shell with the 
fluorescence rate. The refilling rate can be estimated as 
follows: 



re/ 



Here, « is the density of target atoms, cr is the cross- 
section for charge transfer (~ 1CT 14 cm 2 ), and v is the 
relative velocity of the exotic atom with respect to other 
atoms of the medium (< 10 5 cm/s). The typical value 
of the refilling rate is ~ 10 10 s _I for low pressure gases 
and ~ 10 16 s _I for depleted solid and metal 1 18 1. 

Since the Auger transition can take place only if an 
electron occupies a shell state, the time-dependent fill- 
ing condition of the electron in each shell and the re- 
filling rate from outside, including the electron fluores- 
cence transition (de-excitation) from the outer shell to 
the inner shell, Tfl u , needs to be included for a more 
precise calculation in the cascade model with the time 
dependent electron population EH . However, as de- 
scribed below, this will not affect the X-ray yield in the 
low n states since the radiative transition rate dominates 



over the Auger transition rate as n becomes smaller and 
the radiative transition takes place much faster than the 
electron refilling rate. Therefore, we simply estimate 
the modified Auger transition rate, including the elec- 
tron refilling rate and the fluorescence transition rate, 
as: 



YAug,K,mod 
«l,'l-»«2,'2 



ni,h—*n2,h 



T ref 



flu ) 



2.2.2. Radiative Transition 

The radiative transition rate becomes larger than the 
Auger process at a relatively low n state. It can be es- 
timated with a perturbation method and in the dipole 
approximation it follows EOll . 



^Rad 
n\,li—>n2,h 



4e 2 
3/i 4 c 3 



max(/i,/ 2 ) 2 
2/i + l 



AE nun2 = hcRyfiZ 2 



V" 1 



Here, T Rad , , is the radiative transition rate with the 
initial state (n\, l\) and the final state (n.2, h), AE„ um 
is the energy difference between the initial and final 
state, and R v is the Rydberg constant. As seen in the 
equation, the radiative transition rate increases as n de- 
creases (AE nuH2 increases), and becomes the main tran- 
sition process in low n states. The radiative transitions 
dominate for n < 9 for the antiprotonic exotic atom and 
n < 5 for the muonic exotic atom. Note that the radiative 
transitions prefer large An since they are proportional to 
(AE, n „,) 3 , as seen in the equation. However, once the 
cascader reaches the circular state, (n, n - 1), the selec- 
tion rule (Al = ±1) restricts the transition to (n, n - 1) — > 
(n - 1, n - 2). Therefore, we expect a high X-ray yield 
in the low n states, since the cascader is predominantly 
in a circular state at low n. 

2.2.3. Nuclear Capture 

Since the effective Bohr radius for the cascader, ao/fi, 
is much smaller than the Bohr radius, «o, the strong nu- 
clear force interaction between the cascader and the nu- 
cleus can become large in low n states. This may termi- 
nate the de-excitation cascade of the exotic atom before 
it reaches the ground state, since the cascader is cap- 
tured by the nucleus. In particular, the antiproton and 
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the antideuteron annihilate with the nucleus due to the 
nuclear capture and produce pions and protons. The op- 
tical potential between the cascader and the nucleus can 
be estimated as follows Il22l|23l : 



U{r) 



2n 

M* 



1 + \ap(r) 



P(r) 



-(V + iW) 

P(0) 
~ 1 +e'~ 



P(r) 
P(0) 



Here, M* is the reduced mass of the cascader, m# is 
the mass of the nucleon, a is the average complex "ef- 
fective" hadron-nucleon scatting length (experimentally 
determined), and p(r) is the Fermi distribution with the 
parameters p(0) = 0.122 for 3 , c = 1.07 x A 1/3 fm, and 
z = 0.55 fm 11221 |23l |24l. 

The nuclear capture rate can be derived with the per- 
turbation method using the imaginary part of the optical 
potential W, as seen below: 



-■Cap 



-u 



lm(U(r))(R(n u li)) 2 r 2 dr 



2W 

~h~ 



I 



2 ,,2 



(R(nuh) Yr 
l+e l 



dr . 



Here, h is the reduced Planck constant and W is ~ 20 
MeV (experimentally determined J22l |23] EU). Note 
that the energy level of the exotic atom might be slightly 
shifted, due to the strong nuclear force, but the shift is 
small for low and middle Z atoms and negligible com- 
pared with the energy of the atomic X-rays (AE„ U „ 2 ). 

2.3. Parameter Study and Comparison with Experimen- 
tal Data 

A Monte Carlo simulation for the cascade model was 
developed to estimate the X-ray yields of the exotic 
atom. The simulation takes into account all the possi- 
ble Auger transitions including the electron refilling and 
fluorescence transitions, the radiative transitions, and 
the nuclear capture. It starts atn e = 1 (electron K shell), 
and / is determined with the modified statistical distri- 
bution Pi oc (21 + l)e al as discussed above. Cascaders 
are then allowed to cascade until they are captured by 
the nucleus or reach the (1,0) state. The absolute X-ray 
yields, Y, H ^ m , in the low n states (radiative transition 
dominates) were calculated as follows: 



m-i«2-i AT 
/ i= o ij=o v all 



pRad 

r na<- 



Here, the initial and final states are (nuh) and (n2,h) 
(no final state for the nuclear capture), A^ a jj is the num- 
ber of antiprotons simulated in the cascade model and 
N ni j. is the number of antiprotons that cascaded to the 
state (ni,/,-). 

The Monte Carlo simulation was conducted with 
three parameters, a for initial angular momentum dis- 
tribution, r re f for the electron refilling rate, and W for 
the optical potential. Table [T] shows the X-ray yields of 
antiprotonic exotic atoms (Al target) with the different 
values of T re f around the empirical values, T re f = 10 
s~' (a = 0.16 and W = 10 MeV). This indicates, as 
discussed above, the X-ray yields at low n states were 
not affected by the electron refilling rate. The results 
are also consistent with models including the time de- 
pendent electron populatiorQ Table [2] shows the X-ray 
yields of antiprotonic exotic atoms (Al target) with the 
different values of W (a = 0.16 and T ref = 10 16 s" 1 ). 
This also indicates that W affects only one transition, the 
lowest «, as expected. Table [3] shows the X-ray yields 
of antiprotonic exotic atoms (Al target) with the differ- 
ent values of a (W = 10 MeV and T ref = 10 16 s" 1 ). As 
seen in the tables, the X-ray yields are driven mainly 
by a, the initial angular momentum distribution, except 
the last transition is strongly affected by W, the nuclear 
potential. 
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92 keV (5 -> 4) 
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50 keV (6 -> 5) 
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30 keV (7 -> 6) 
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Table 1 : X-ray yields of the antiprotonic exotic atom (Al target) with 
different values of T. 
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50 
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92 keV (5 -> 4) 


89% 


67% 


46% 


35% 


22% 


50 keV (6 -> 5) 


82% 


82% 


82% 


81% 


82% 


30 keV (7 -> 6) 


69% 


69% 


69% 


68% 


69% 



Table 2: X-ray yields of the antiprotonic exotic atom (Al target) with 
different values of W. 



Additionally, our cascade model was compared with 
the experimental data for muonic exotic atoms, which 
are widely measured in experiments ll25ll . The nuclear 
absorptions are not seen in the muonic exotic atoms 
(except for high Z targets) and therefore, there is only 
one parameter, a, to control the X-ray yields at low n 



1 private discussion with Dr. Takahisa Koike (RIKEN Japan) 
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a 





0.08 


0.16 


0.24 


92 keV (5 -> 4) 


41% 


58% 


67% 


71% 


50 keV (6 -> 5) 


46% 


69% 


82% 


88% 


30 keV (7 -> 6) 


37% 


56% 


69% 


75% 



Table 3: X-ray yields of the antiprotonic exotic atom (Al target) with 
different values of a. 



states. Table |4] shows the comparison with the exper- 
imental data and a cascade model developed by Vogel 
and Hartmann for the muonic exotic atoms [25, 26, 27]. 
The parameters used here were W = MeV (no nuclear 
absorption), T ref = 10 16 s"\ and a = 0.18, -0.16, -0.01 
for Al, Fe and Au targets. 



Transition 


exp 


our model 


model in [26, |27l 


Al (2 -> 1) 


80% 


80% 


80% 


(3-»2) 


63% 


63% 


60% 


(4-»3) 


34% 


41% 


42% 


Fe (2 -> 1) 


72% 


71% 


74% 


(3-»2) 


44% 


49% 


45% 


(4-»3) 


33% 


35% 


33% 


Au (2 -> 1) 


90% 


94% 


95% 


(3-»2) 


80% 


84% 


84% 


(4-»3) 


76% 


73% 


76% 



Table 4: Experimental data and cascade models for X-ray yields of 
the muonic exotic atoms with Al, Fe and Au targets 



As seen above, X-ray yields at low n states in our cas- 
cade model are in good agreement with both experimen- 
tal data (muonic exotic atoms) and other cascade models 
(both muonic and antiprotonic exotic atoms) with a time 
dependent electron population. 

3. Accelerator Test at KEK 

3.1. Overview of Accelerator Test 

The KEK facility is located north of Tokyo, in 
Tsukuba, Japan. During the course of the experiments 
the proton synchrotron produced an 8 GeV (up to 12 
GeV) proton beam in the main ring. The H ion source 
generated in the plasma chamber was injected into the 
pre-injector, followed by the linac, booster synchrotron 
and main ring and accelerated to 750 keV, 40 MeV, 500 
MeV and 8 GeV, respectively. Our experiment was per- 
formed at the 7r2 secondary beam line, which delivers 
copious particles including antiprotons generated by the 
proton beam hitting an internal target in the main ring 
(see Figure p). 

The beam test was conducted at KEK in 2004 and 
2005 to verify the GAPS original concept described in 




Figure 7: KEK beam line 

(2) and measure the X-ray yields of the antiprotonic ex- 
otic atom with several different target materials. The 
results constrained the parameters in the cascade model 
described in Section |2] This also allowed us to extend 
the cascade model to any exotic atoms and estimate the 
X-ray yield of the antiprotonic and antideuteronic exotic 
atoms in the GAPS experiment. 



3.2. Experimental Setup 




Figure 8: KEK experimental setup 

The GAPS instrument in the KEK test was composed 
of a TOF, degraders (lead brick and sheets), shower 
counters, a target and X-ray detectors. The antipro- 
tons in the beam were first identified by the TOF sys- 
tem, since antiprotons are slower than the other parti- 
cles in the beam. The degrader slowed down antipro- 



8 



tons and stopped them in the target material where they 
formed an excited antiprotonic exotic atom. Atomic X- 
rays and charged particles are emitted in the decay of 
the exotic atom as discussed in the previous chapter. A 
Sodium Iodide doped with Thallium, Nal(Tl), detector 
array was installed around the target material and de- 
tected the atomic X-rays and pions. The shower coun- 
ters monitored the energy deposited by the particles in 
the beam and distinguished antiprotons from other par- 
ticles, including the in-flight annihilation products. 

While gaseous targets and a few liquid targets were 
used in 2004, liquid and solid targets were tested in 
2005, since they are simpler to implement in the realis- 
tic design for the balloon experiment. The actual picture 
and schematics of the experimental setup in 2005 are 
shown in Figures [8] and [9] Figure 10 shows an unfolded 
view of the cylindrical detector module and a typical 
stopped antiproton event for the Al target. Numbers are 
the energy of the stopped X-rays and n* indicates a pion 
hit. 



Nal Energy [keV] 



Dl'lji'.ilLt 




Figure 9: Experimental setup at KEK in 2005. 



3.2.1. Beam Profile 

The momentum of the beam was controlled and fo- 
cused by dipole and quadrupole magnets, while the 
momentum spread was controlled by a shutter. The 
particles were delivered in 1.5 s long spills, and each 
spill was separated by a 4 s interval. A momentum 
of 1 GeV/c was used in all of the 2005 measurements. 
Up to several GeV/c, the antiproton flux from the n2 
beam line increases with increasing momentum; how- 
ever, losses due to annihilation in the degrader increase 
for the thicker degraders required to stop higher mo- 
mentum antiprotons. Around 1 GeV/c was found to 
provide the highest rate of antiproton stops in our tar- 
get during our 2004 measurements. The beam spill with 
a momentum of 1 GeV/c contained about 20-30 antipro- 
tons, 10 5 7r~, and a somewhat smaller number of K~ and 
e~, as measured in the 2004 experiment, and these num- 
bers were consistent with the data sheets provided by 



Time |ns] 
E[MeV| 



Figure 1 0: A typical antiproton event for the Al target. Numbers are 
the energy of the stopped X-rays and x* indicates a pion hit. 



KEK. 

The spatial beam profile at the P0 counter was mea- 
sured by changing the last dipole magnet, which con- 
trolled the horizontal direction, and the height of the re- 
mote controlled table (for the vertical direction). The 
measured beam profiles at P0, PI and P2 were used as 
input in the GEANT4 simulationFhogether with a TUR- 
TLE beam line optics ray trace rjto simulate the beam 
profile, divergence and momentum bite. 

3.2.2. Time of Flight System 

The TOF system in the KEK 2005 test was composed 
of 6 scintillation counters, P0-P5. The TOF timing, the 
travel time of the incoming particle between the P0 and 
Pi (i = 1, 2, 3, 4, 5) counters, allowed us to identify 
the incoming particle, since all the particles in the beam 
had a fixed momentum and the antiprotons were much 
slower than the other lighter particles (see below). The 
P0, P2, P3 and P4 counters had a dimension of 12 cm x 
12 cm and a thickness of 1.0 cm, while the PI and P5 
counters had a thickness of 0.2 cm. The paddles were 
coupled to the light guide and then to the 2 inch fast 
photomultiplier tube (Photonis XP2020). A high volt- 
age of ~ -1800V was applied to the PMT base (Photonis 
S5632). 

The P0, PI, and P2 counters were used for tim- 
ing only, while the P3, P4 and P5 counters were used 



2 GEometry ANd Tracking, a toolkit for the simulation of the pas- 
sage of particles through matter, developed by CERN. 

3 PSI Graphic Turtle Framework by U. Rohrer based on a CERN- 
SLAC-FERMILAB version by K.L. Brown et al. 
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for both timing and energy deposition. The timing at 
each counter was measured relative to the accelerator 
beam structure by passing the signal from the last diode 
through a fast timing preamplifier (Ortec VT120b), fol- 
lowed by a constant fraction discriminator. The time 
of flight (TOF) between the PO counter and the P1-P5 
counters were measured using time- to-analog convert- 
ers (TAC, Canberra 2020). The dE/dX energy deposit 
was characterized using the signal from the PMT an- 
ode passed through a preamplifier (Camberra 2005) fol- 
lowed by a spectroscopy amplifier (Ortec 452). 

3.2.3. Shower Counter 

The four shower counters, S1-S4, were installed be- 
hind the main degrader in 2005 (see Figure[9]i, and each 
of them had a dimension of 12 cm x 12 cm x 0.5 cm. 
A 0.25 inch lead sheet was sandwiched between every 
pair of counters to slow down the incoming particle. 
Each counter was coupled to the light guide and PMT 
(Photonis 2042 and 2072 PMT). The shower counter 
allowed us to distinguish antiprotons from other parti- 
cles by measuring the dE/dX energy loss, since non- 
relativistic slow antiprotons deposit more energy than 
relativistic particles such as n~ . 

The veto counters (6 cm wide, 1 mm thick ribbon 
scintillation fibers, coupled to a Hamamatsu R1942A 1 
inch PMT) were installed between the target and the 
X-ray detectors. They were designed to monitor the 
off-axis antiprotons hitting the detector and the frame 
without stopping in the target material. However, since 
the energy resolution of these counters was relatively 
coarse, it was difficult to uniquely identify potential off- 
axis antiproton interactions from annihilation products 
produced in the target. 

3.2.4. X-ray Detector 

The X-ray detectors were 128 Nal(TI) crystals (1 inch 
x 1 inch x 5 mm). The Nal(TI) detector emits scintil- 
lation light proportional to the deposited energy, ~ 40 
photons/keV. Since the Nal(Tl) is a relatively high Z 
material, up to 300 keV X-rays can be photo-absorbed 
in the 5 mm thick crystal. Each crystal is coupled to a 
Hamamatsu 1 inch PMT (R 1924 A) on the back surface. 
The wavelength of the scintillation light is ~ 410 nm, 
where the quantum efficiency of the PMT has a peak. 
Every 8 crystals and PMTs, separated from each other 
by 1.5 inch, are mounted in a tightly sealed steel housing 
with a 0.125 mm Al window. Each PMT is connected 
to the custom made PMT base and ~ -800V HV was ap- 
plied. The preamplifier was mounted inside the housing 
and the gain for each detector was controlled externally. 



Sixteen sets of detectors were mounted around the tar- 
get as seen in Figure 1 1 




Figure 1 1 : Picture for KEK X-ray array. 

3. 2. 5. Target Material 

In 2005, four target materials were chosen based on 
the energy of the atomic X-rays in their antiprotonic 
exotic atom, which needed to be in the useful energy 
range of the X-ray detector, 25 keV < E < 300 keV. 
The detectable antiprotonic atomic X-rays for each tar- 
get tested in KEK are shown in Table [5] 



Table 5: Antiprotonic atomic X-rays for each target (25 keV < E < 
300 keV) 



Target 


X, 


X 2 


x 3 


X 4 


x 5 


x 6 


Al 


92keV 


50keV 


30keV 








S 


139 keV 


76keV 


46keV 


30keV 






CI 


86keV 


52keV 


34keV 


23 keV 






Br 


145 keV 


99keV 


71 keV 


52keV 


41 keV 


31 keV 



Al (Aluminum wool), S (Sulfur), CBr 4 (Tetrabro- 
momethane) and CCI4 (Carbon tetrachloride) targets 
were tested in 2005 and we will focus on the Al (Z = 
13) and S (Z = 16) targets in this paper to estimate the 
X-ray yields for the Si (Z = 14) target for the GAPS bal- 
loon experiment. It is also more challenging to analyze 
the data for the CBr4 and CCI4 targets since they are 
compounds and many atomic X-rays can be produced 
in the small energy region. The Al wool was filled into 
two 1 mm thick plastic bottles, each with a diameter of 
12 cm and 22 cm in length, and the average density was 
~ 0.111g/cm 3 . The target holder for the Sulfur powder 
was framed with Al pipes of diameter 12 cm cut at a 
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45 degree angle, and both openings were covered with 
1 mm thick plastic sheets (see Figure 12 1. The hold- 
ers were placed onto two guided rails to minimize the 
blockage of X-rays from the target. 





9 10 11 12 

degrader thickness [cm] 



Figure 12: Sulfur target geometry 



Figure 13: Number of hits at each counter vs. degrader thickness 
(cm). The counts were normalized to the counts at P0 counter. 



3.2.6. Degrader Thickness and Range Curve 

Since antiprotons in the beam were too energetic to 
stop in the target, a combination of active and passive 
degraders were used to slow down the antiprotons be- 
fore they entered the GAPS target region (see Figure|9]). 
The optimized total thickness of degrader was estimated 
by measuring the number of events at the P4 counter 
(just before the target) with different thicknesses of de- 
grader. Since the number of antiprotons in the beam 
was very small, in order to have better statistics, we 
used positively charged beam (protons and n + ) with the 
same magnet settings for the beam except for the po- 
larity. Figure [13] shows the number of protons at each 
counter, normalized with the number of protons at the 
P2 counter. The GEANT4 simulation result at the P4 
counter was also shown in the figure, taking into ac- 
count the uncertainty on the thickness and density for 
each lead brick and sheet (~ 3%). The data and the 
simulation result are in good agreement with each other 
and the number of protons at the P4 counter rapidly de- 
creased as the thickness of the degrader increased to ~ 
10.3 cm. This implies that there were many slow pro- 
tons present at the P4 counter with this thickness and 
thus we decided to use a total thickness of the degrader 
of 10.3 cm in the experiment. 
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Figure 14: TOF timing at TAC1, TAC2, TAC3, TAC4, and TAC1 vs. 
TAC2. 



3.3. Antiproton Selection 

Since the momentum of the beam was set as 1 GeV/c 
by the dipole and quadrupole magnets, the antiprotons 
in the beam can be distinguished from other particles 
by the velocity, The TOF timing and energy deposit 
in the plastic scintillator for antiprotons are larger than 
other particles in the beam (mainly pions) since (3 for 



antiprotons is smaller and dE/dx energy loss is propor- 
tional to ~ j£T 2 . 



Figure [T4| shows the TOF timing at TAC1, TAC2, 
TAC3 and TAC4 (between the P0 and Pi (i = 1, 2, 3, 4) 
counters). A plot for TAC1 vs. TAC2 was also shown. 
Two peaks, relativistic pions (pre-scaled) and antipro- 
tons, are seen at each plot. The selected cuts, peak ± 1 
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ns, were applied for each TOF timing to select antipro- 
ton events. Since the P5 counter was placed ~ 60 cm 
away from the P4 counter, the cut on the TAC5 was set 
as "(TAC4 lower limit + 2 ns) or no hit". Tables[6]shows 
the applied cuts for each TOF timing. 





lower limit 


upper limit 


TAC1 


14.0 ns 


16.0 ns 


TAC2 


29.0 ns 


31.0 ns 


TAC3 


30.5 ns 


32.5 ns 


TAC4 


30.5 ns 


32.5 ns 


TAC5 


> 32.5 ns or no hits 



Table 6: Antiproton selection cuts on each TOF timing 




Figure 15: dE/dX energy deposit in the SI, S2, S3 and S4 counters. 



Figure 15 shows the dE/dX energy deposit in the 
shower counters. Energy was calibrated with the rela- 
tivistic pions in the beam and the GEANT4 simulation 
with the actual beam profile. Two peaks, relativistic pi- 
ons (pre-scaled) and antiprotons, are also seen at each 
plot. The antiproton selection cuts were applied for each 
dE/dX energy deposit. The applied cuts for each dE/dX 
energy deposit are shown in Table [7] Note that we were 
not able to set the upper limits on the cuts for the P3 
and P4 dE/dX energy deposits since the signals were 



saturated at E > 10 MeV. Figures 16 and 17 show the 
spectrum in the detector before and after applying the 
cuts. The X-ray lines are seen in the spectrum with the 
cuts. 

3.4. Background Model 



As seen in Figure 17 the cuts applied to the TOF 
timing and dE/dX energy deposit provide excellent an- 
tiproton selection in the original beam. Thus, the main 





lower limit 


upper limit 


SI 


1.8 MeV 


3.2 MeV 


S2 


2.2 MeV 


4.2 MeV 


S3 


2.2 MeV 


4.2 MeV 


S4 


2.6 MeV 


5.0 MeV 


P3 


8.0 MeV 




P4 


8.0 MeV 





Table 7: Antiproton selection cuts on each dE/dX energy deposit 




250 
energy [KeV] 



Figure 16: Energy spectrum in the X-ray detector without applying 
any cuts. 




Figure 17: Energy spectrum in the X-ray detector with the cuts shown 
in Table|7] 



background is due to the annihilation products of the 
stopped antiprotons, which can develop an electromag- 
netic shower in the target and the detector frame. Most 
of the antiprotons in the beam were annihilated in the 
degrader and the annihilation products developed the 
electromagnetic shower in the target and the detector 
frame. Therefore, the background model was obtained 
from the experimental data with cuts on the TOF timing 
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at TAC1 (between the PO and PI counters) and TAC2 
(between the PO and P2 counters) to evaluate the elec- 
tromagnetic shower generated by the annihilation prod- 
ucts. Additionally, another background model was ob- 
tained from a GEANT4 simulation with the stopped an- 
tiprotons around the target. Figures 18 and 19 show 



the spectrum for each background model. Both mod- 
els show the similar shape of the background and there 
is a wide peak around 100 keV. The former background 
model (called KEK BG) uses the actual data and detec- 
tor geometry, however, it is for the antiproton annihi- 
lation in flight in the degrader. On the other hand, the 
latter background model (called GEANT4 BG) was ob- 
tained from the GEANT4 simulation with the stopped 
antiprotons around the detector, but we do not have a 
good way to estimate the systematic uncertainties asso- 
ciated with the physics model and the detector geom- 
etry. Therefore, in the following sections, we will use 
both models to fit the data to obtain the X-ray yields of 
the exotic atom. 




250 

energy [KeV] 



Figure 18: The background model for the Al target, obtained from the 
experimental data with the cuts on the TOF timing at TAC 1 and TAC2 
applied. 



3.5. Atomic X-ray Spectrum 

Since antiprotons were not able to be tracked after 
hitting the P4 counter, GEANT4 simulation was used 
to predict the stopped position of the incoming antipro- 
tons. The simulation was also used to estimate the en- 
ergy spectrum in the detector for each atomic X-ray, tak- 
ing into account all the X-ray interactions before reach- 



ing the detector. Figure 20 shows the expected energy 
spectra in the detector for each atomic X-ray with the 
Al target. It is normalized to the counts per exotic atom 
with 100% X-ray yield, and the green, blue, and red 
lines represent 30 keV, 50 keV, and 92 keV X-rays. The 
dashed lines are the simulation results without the de- 
tector response and the solid lines are the spectra with 
the detector response (7% FWHM at 1 MeV). Figure |2T 
is the same for the S target (30 keV for purple, 46 keV 
for green, 76 keV for blue and 139 keV for red lines). 
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Figure 20: Expected energy spectra in the detector for each atomic 
X-ray with the Al target. The green, blue, and red lines represent 
simulation results for 30 keV, 50 keV, and 92 keV lines, and the solid 
lines are the spectra with the detector response. It is normalized to the 
counts per exotic atom with 100% X-ray yield. 




energy deposit [KeV] 



Figure 19: The background model for the Al target, obtained from the 
GEANT4 simulation with the stopped antiprotons around the target. 



3.6. X-ray Yields of Antiprotonic Exotic Atom 

The absolute X-ray yields (probability to emit an 
atomic X-ray per exotic atom) for antiprotonic exotic 
atoms, Y, were estimated by fitting the data with the 
background model and the expected energy spectra for 
each atomic X-ray in the detector as below. 

- a B c ■ /bc + ^ <k ■ fi 

Yi = at/"*** 

Here, a's are the fit parameters, i is for the z'-th atomic 
X-ray, /'s are the spectra in the detector as discussed 
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Figure 2 1 : Expected energy spectra in the detector for each atomic X- 
ray with the S target. The purple, green, blue, and red lines represent 
simulation results for 30 keV, 46 keV, 76 keV and 139 keV lines, 
and the solid lines are the spectra with the detector response. It is 
normalized to the counts per exotic atom with 100% X-ray yield. 




Experimental data 

Al X-rays (30keV. 50keV. 92keV) 

KEK BG 

Fitting result 



250 

energy [keV] 



Figure 23: The data for the Al target fitted with the background model 
obtained from the experimental data (blue) and the expected X-ray 
spectrum for each antipro tonic X-ray (red). 



above, and Np sto? is the number of stopped antiprotons 
in the target. The parameter a,- also denotes the number 
of the atomic X-rays emitted from the exotic atom. 

3.6.1. Al Target 




Experimental data 

Al X-rays (30keV, SOkeV. 92keV) 

GEANT BG 

Fitting result 



\ 



250 
energy [keV] 



Figure 22: The data for the Al target fitted with the GEANT4 back- 
ground model (blue) and the expected X-ray spectrum for each an- 
tiprotonic X-ray (red). 



Figures 22 and 23 show the fitting results with the 
two different background models for the Al target. The 
solid black, blue and red lines represent the experimen- 
tal data, the background model and the three atomic X- 
rays respectively and the green solid line is the sum of 
the background model and the X-rays. Table [8] shows 
the X-ray yields for each atomic X-ray including the fit- 
ting error (Y + AY). High absolute yields ~ 80% were 
seen for all three transitions in both background models. 
The nuclear absorptions were not seen in the n — 5 — > 4 



p-Al 


KEK BG 


GEANT4 BG 


92 keV (5 -> 4) 


90% ±11% 


81% ± 10% 


50 keV (6 -> 5) 


76% ± 8% 


85% ± 8% 


30 keV (7 -> 6) 


84% ±11% 


90% ±11% 


reduced-^ 2 


1.07 


1.09 



Table ! 
els. 



: The fitting results with the error for the two background mod- 



transition, but were seen in the n — 4 — > 3 transition. 
The systematic uncertainty due to the detector response 
(AFWHM ~ ±1% at 1 MeV) and the offset of the energy 
calibration (± 1 keV) was estimated as ~ 7%, while the 
uncertainty due to the background model is ~ 4%. As a 
result, the absolute yields for the atomic X-rays (Y+AY) 
were obtained as 87% ± 13% for 30 keV, 81% ± 11% 
for 50 keV, 86% ± 13% for 92 keV (see Table [TO}. 

3.6.2. S Target 

Since some of the antiprotons may stop in the Al win- 
dow/frame around the target, seven atomic X-rays (three 
from the exotic atoms with the Al window/frame and 
four from the S target) can be produced in the small en- 
ergy region. Therefore, in order to simplify the fitting, 
we constrained the three atomic X-rays for the S target, 
30 keV (n = 8 -> 7), 46 keV (n = 7 -> 6) and 76 keV 
(n — 6 — > 5), to have the same absolute yields. This is 
consistent with theoretical expectation for the S target 
and as predicted and measured for the Al target. Ad- 
ditionally, the X-ray yields for the antiprotons stopped 
in the Al window/frame are constrained with the value 
obtained for the Al target. 

Figures [24] and [25] show the fitted results with the two 
different background models for the S target. As seen in 
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Experimental data 

S X-rays (30keV. 46keV, 76keV. 139keV) 
GEANT BG 
Fitting result 




Figure 24: The data for the S target fitted with the GEANT4 back- 
ground model (blue) and the expected X-ray spectrum for each an- 
tiprotonic X-ray (red ). 




Experimental data 

S X-rays (30keV, 46keV, 76keV, 139keV) 
KEK BG 
Fitting result 



50 100 



150 200 



250 
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Figure 25: The data for the S target fitted with the background model 
obtained from the experimental data (blue) and the expected X-ray 
spectrum for each antiprotonic X-ray (red). 



p-S 


KEK BG 


GEANT4 BG 


139 keV (5 -> 4) 


59% ± 19% 


20% + 20% 


76 keV (6 -> 5) 


72% ± 17% 


66% ± 16% 


46 keV (7 -> 6) 


72% ± 17% 


66% + 16% 


30 keV (8 -> 7) 


72% ± 17% 


66% + 16% 


reduced-^ 2 


1.00 


1.39 



Table 9: The fitting results with the error for the two background mod- 
els. 



3. 7. Comparison with Cascade Model 

As discussed in Section [2] the cascade model has 
been developed to estimate the X-ray yields of the ex- 
otic atoms at low n states. The yields are mainly de- 
termined by the parameter a, initial angular momentum 
distribution, while the last transition rate is strongly de- 
pending on W, nuclear potential (see Section [2j. The 
cascade model with the parameters, a — 0. 16 and W = 5 



MeV (T, 



ref 



10 s ) is quite consistent with the ex- 



perimental data for the Al target as seen in Table 10 The 
absolute yield for the 30 keV X-ray is higher than we 
expected in the cascade model, but this could be due to 
the atomic X-rays produced from other materials around 



the detector. Table 1 1 shows the X-ray yields of the ex- 
perimental data and the cascade model for the S target, 
which is also in good agreement with the experimental 
data. Parameters used here are, a = 0.16 and W = 5 
MeV <T re/ = 10 16 s- 1 ). 



p-Al 


Experiment 


Cascade Model 


92 keV (5 -> 


4) 


86% ± 13% 


76% 


50 keV (6 -> 


5) 


81% ± 11% 


82% 


30 keV (7 -> 


6) 


87% + 13% 


69% 



Table 10: X-ray yields and experimental data for antiprotonic exotic 
atom with Al target. 



Table [9j the reduced chi-square for the GEANT4 back- 
ground model became much higher than the one for the 
KEK model. This is due to the difficulty in modeling the 
complex S target geometry in the GEANT4 simulation. 
Thus, the X-ray yields of the antiprotonic exotic atom 
with the S target were estimated with the KEK back- 
ground model as follows, by taking into account all the 
systematic uncertainties discussed above; 72% + 18% 
for 30 keV, 72% ± 18% for 46 keV, 72% ± 18% for 76 
keV and 59% ± 21% for 139 keV (see Table [l]}. High 
absolute yields were also seen in all the transitions ex- 
cept for the n — 5 — > 4 transition, which is due to the 
nuclear absorption and consistent with the result shown 
in 11281. 



p-S 


Experiment 


Cascade Model 


139 keV (5 -> 4) 


59% + 21% 


48% 


76 keV (6 -> 5) 


72% ±18% 


78% 


46 keV (7 -> 6) 


72% + 18% 


70% 


30 keV (8 -> 7) 


72% + 18% 


54% 



Table 1 1 : X-ray yields and experimental data for antiprotonic exotic 
atom with S target. 



3.8. Prediction for Antiprotonic and Antideuteronic Ex- 
otic Atom for Si Target 

The cascade model was extended to the antiprotonic 
and antideuteronic exotic atom for a Si target and other 
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materials in the GAPS instrument to estimate the an- 
tideuteron sensitivity. Since parameters are not strongly 
correlated with the atomic number, the same parameters 
are used for the Si (Z = 14) target as used for the Al (Z 
= 13) and S (Z = 16) targets. Table 12 shows the re- 
sult for the antiprotonic exotic atom with the Si target 
(a = 0.16, W = 5MeV,r re/ 



10 16 s- 1 ). 



p-Si 


Cascade Model 


Estimation in 1 2 1 


106 keV (5 -> 4) 


70% 


50% 


58 keV (6 -> 5) 


81% 


50% 


35 keV (7 -> 6) 


70% 


50% 



Table 12: Cascade model for X-ray yields of the antiprotonic exotic 
atom with the Si target. 

The X-ray yields for the antideuteronic exotic atom 
with a Si target were also estimated by simply changing 
the optical potential, Wj ~ 2Wp = 10 MeV, as shown 
in Table [B] (a = 0.16, T ref = 10 16 s" 1 ). It was also 
estimated for higher values of W, however, the result 
does not affect the GAPS antideuteron sensitivity since 
the nuclear capture only takes place at n — 5 or n — 6 
for extremely high values of W, and the corresponding 
atomic X-ray energies (206 keV, 112 keV) are too high 
to be detected in the GAPS detector. The result indicates 
an increase of GAPS sensitivity to the antideuterons [0 
since the X-ray yields for the GAPS experiment were 
previously assumed to be ~ 50 % in |H. 



the antiprotonic and muonic exotic atoms. The model 
allows us to estimate the X-ray yields of the antipro- 
tonic and antideuteronic exotic atoms formed with any 
materials in the GAPS instrument and the X-ray yields 
for antiprotonic and antideuteronic exotic atoms with a 
Si target were estimated as ~ 80%. This is higher than 
previously assumed in [2], indicating the increase of the 
GAPS antideuteron sensitivity. Consequently, we ob- 
tained the GAPS antideuteron sensitivity as 1.3 xl0~ 6 
[rrrV'sr-'CGeV/n)- 1 ] for a 105 day, 3 flight balloon 
program 0, which has a strong potential to detect an- 
tideuterons produced by dark matter. 
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d-Si 


Cascade Model 


Estimation in t2l 


206 keV (5 -> 4) 


24% 




112keV (6-> 5) 


83% 




67 keV (7 -> 6) 


82% 


50% 


44 keV (8 -> 7) 


76% 


50% 


30 keV (9 -> 8) 


65% 


50% 
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atom with Si target. 
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